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ABSTRACT 

Bacterial Hfq is a protein that plays an important 
role in the regulation of genes in cooperation with 
sRNAs. Escherichia coli Hfq (EcHfq) has two or more 
sites that bind RNA(s) including U-rich and/or the 
poly(A) tail of mRNA. However, functional and struc- 
tural information about Bacillus subtilis Hfq (BsHfq) 
including the RNA sequences that specifically 
bind to it remain unknown. Here, we describe 
RNA aptamers including fragment (AG) 3 A that are 
recognized by BsHfq and crystal structures of 
the BsHfq-(AG) 3 A complex at 2.2 A resolution. 
Mutational and structural studies revealed that the 
RNA fragment binds to the distal site, one of the two 
binding sites on Hfq, and identified amino acid 
residues that are critical for sequence-specific 
interactions between BsHfq and (AG) 3 A. In particu- 
lar, R32 appears to interact with G bases in (AG) 3 A. 
Poly(A) also binds to the distal site of EcHfq, but the 
overall RNA structure and protein-RNA interaction 
patterns engaged in the R32 residues of BsHfq- 
(AG) 3 A differ from those of EcHfq-poly(A). These 
findings provide novel insight into how the Hfq 
homologue recognizes RNA. 

INTRODUCTION 

Gene expression is controlled in bacteria through re- 
sponses to environmental changes that affect growth. 
Hfq is an important component in the regulation of 
gene expression in cooperation with sRNAs. Hfq (also 
known as HF-I) was originally identified as an 
Escherichia coli host factor required for the RNA 



replication of phage Q(3 (1) and further studies revealed 
that Hfq in Gram-negative bacteria functions as a 
post-transcriptional regulator by interacting with various 
RNAs. Escherichia coli Hfq (EcHfq), which is the most 
studied among this family of proteins, promotes inter- 
actions between small untranslated RNA regulatory mol- 
ecules (such as OxyS, DsrA, Spot42, RyhB and SgrS 
sRNAs) and their target mRNAs (2-6). Escherichia coli 
Hfq interacts with many target mRNAs as an RNA chap- 
erone (7,8) and stimulates polyadenylation that is 
catalyzed by poly(A) polymerase I, by binding to the 
poly(A) tail (9-11). Furthermore, EcHfq is involved in 
RNA stability control by protecting RNA from degrad- 
ation by the 3'-5' exoribonucleases, RNase II and 
PNPase (12,13). 

About half of all sequenced Gram-negative (such as 
E. coli spp., Salmonella enterica, Pseudomonas aeruginosa, 
Neisseria meningitidis and Vibrio cholerae) and 
Gram-positive (e.g. Staphylococcus aureus, Listeria 
monocytogenes and Bacillus subtilis) bacteria express Hfq 
[reviewed in Ref. (14)]. Growth rates of Hfq deletion 
strains of E. coli K-12 and P. aeruginosa Ol are reduced 
(15,16). In addition, growth of the Hfq deletion strain of 
TV. meningitidis is impaired in nutrient-rich media and 
proteomics analysis has revealed that the expression of 
28 genes is affected in this strain (17). Moreover, Hfq is 
required for intestinal colonization by V. cholerae (18). On 
the other hand, little is known about the function of Hfq 
in Gram-positive bacteria. The growth of S. aureus and 
L. monocytogenes Hfq deletion strains is not affected 
(19,20). Although Hfq in S. aureus (SaHfq) does 
not play important role in the stress response, RNA sta- 
bility or exoprotein expression (19), L. monocytogenes 
Hfq controls the expression of numerous stress- and 
virulence-associated genes and binds to sRNAs (20,21). 
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The deletion of B. subtilis Hfq (BsHfq), which is 
encoded by the ymaH gene, does not affect either the 
growth or sporulation (22). BsHfq binds to SRI sRNA 
and to ahrC mRNA and the interaction of SRI sRNA 
with its primary target ahrC mRNA causes the inhibition 
of ahrC mRNA translation. However, BsHfq neither sta- 
bilizes SRI sRNA, nor promotes complex formation 
between SRI sRNA and ahrC mRNA (23,24). 
Alternatively, BsHfq is required to activate ahrC mRNA 
translation (24). 

Crystal structures for several bacterial Hfqs and 
archaeal Hfq-like proteins show that Hfq forms a 
homohexameric ring consisting of Sm-like folds and two 
RNA-binding sites at proximal/distal locations in both 
sides of the ring (25-30). A co-crystal structure of SaHfq 
in complex with AU 5 G has revealed that the RNA is 
recognized by residues at its proximal site (two loops 
between (32-(33 and (34-P5), although evidence that RNA 
binds to the distal site of SaHfq has not been uncovered 
(25). Site-directed mutagenesis of EcHfq has revealed that 
the two binding sites recognize their respective RNAs, 
namely poly(A) and U-rich RNA; the poly(A) tail of the 
target mRNA interacts with residues (Y25 and 130) on 
the distal site of EcHfq and short U-rich RNA (AU 5 G) 
recognizes residues (Y55, K56 and H57) on the proximal 
site (31). Moreover, a co-crystal structure of EcHfq in 
complex with poly(A) RNA has recently been described. 
The structure shows that the poly(A) RNA binds to the 
distal site (on loop between (31-(32, (34 and (35) and the 
poly(A)-binding site comprises an adenosine specificity 
site, a purine nucleotide selectivity site, and a sequence- 
non-discriminating RNA entrance/exit site (30). 

In this article, we present the RNA sequence motif that 
is recognized by BsHfq and the crystal structure of 
BsHfq-RNA complex. We first performed a systematic 
evolution of ligands by exponential enrichment (SELEX) 
experiment and identified a single stranded AG repeat 
sequence in specific RNA sequences that are recognized 
by BsHfq. Structures of BsHfq in complex with the RNA 
fragment (AG) 3 A were determined at a resolution of 2.2 A 
by X-ray crystallography using molecular replacement. 
The overall structure of BsHfq is a hexameric ring 
comprising a Sm-like motif, and the RNA fragment 
(AG) 3 A is bound to the distal site of BsHfq. This 
binding mode was confirmed by gel electrophoresis 
mobility shift assays (EMSA) with site-directed mutagen- 
esis. The present results provide novel insight into the rec- 
ognition pattern of RNA with AG repeats located at the 
distal site of Hfq. 

MATERIALS AND METHODS 

Plasmid construction 

We constructed pQE60Hfq for C-terminal 6x His fusion 
BsHfq (BsHfq-His) overexpression by inserting a 
PCR-amplified ymaH gene into the Ncol-Bglll digested 
pQE60 expression vector (Qiagen). The PCR product was 
obtained using ymaH-Ncol and ymaH-Bglll primers 
(Supplementary Table SI). Site-directed mutagenesis was 
performed using a QuikChange site-directed mutagenesis 



kit (Stratagene). Single specific mutations were introduced 
into BsHfq-His protein using pQE60Hfq as a template 
DNA together with mutagenic primers (see primer list in 
Supplementary Table SI). All constructs were verified by 
DNA sequencing. 

Protein overexpression and purification 

The expression of BsHfq-His in E. coli M15/pREP4 was 
induced for 5 h in the presence of 2 mM IPTG. The cells 
were harvested by centrifugation at 4000 g at 4°C for 
15min. Wet cells were suspended in lysis buffer (10 mM 
Tris-HCl, 100 mM NaH 2 P0 4 , 8M Urea, pH 8.0), lysed 
during incubation at room temperature for 30min and 
centrifuged at 15 000g at room temperature for 15min. 
The supernatant was incubated with His-select™ Nickel 
Affinity gel (Sigma) at room temperature for 1 h, and 
eluted with elution buffer A (10 mM Tris-HCl, lOOmM 
NaH 2 P0 4 , 8M Urea, 200 mM Imidazole, pH 6.3). The 
purity of the BsHfq-His protein was analyzed by 15% 
SDS-PAGE. Eluted samples were dialyzed against 
20 mM Tris-HCl buffer (pH 7.5) containing 500 mM 
NaCl and 50% glycerol. The present study assesses 
hexameric conformations of purified single specific 
BsHfq mutations by gel filtration chromatography on 
Superdex™ 200 10/300 GL columns (GE Healthcare) 
equilibrated with binding buffer A (10 mM Tris-HCl, 
pH 7.5, 50 mM NaCl, 50 mM KC1, 1 mM MgCl 2 ). 
Mutant BsHfqs (48 jig) with lysozyme (50 jig) were 
eluted with the same buffer. We purified GST-tagged 
BsHfq and removed the GST-tag using PreScission 
protease for crystal structure analysis and EMSA experi- 
ments with short RNAs. The purified BsHfq contained an 
extra GPLGS sequence at the N-terminus as described 
previously (32). 

SELEX experiment 

RNA selection experiments proceeded as described (33) 
using the template 5 ; - TAATACGACTCACTATA GGG 
ACACAATGGACG - N30 - TAACGGCCGACATGA 
GAG-3' where N30 represents 30 random nucleotide pos- 
itions (T7 RNA polymerase promoter is underlined) and 
the primer sequence, 5 / -CTCTCATGTCGGCCGTTA-3 / . 
We enzymatically synthesized RNA pools using 
AmpliScribe™ T7 High Yield Transcription kit 
(Epicentre Biotechnologies). After the transcribed RNAs 
were incubated with DNase I at 37°C for 15min, RNAs 
were separated by 6% PAGE under denaturing condition 
with 8 M urea and then purified from the gel. 

Aptamers that bind to BsHfq-His were obtained after 
nine rounds of selection. The first six rounds comprised 
filtering the Hfq-RNA ap tamer complex through 0.22 jim 
nitrocellulose filter (Millipore). The samples were mixed in 
binding buffer A. The filters were washed with the binding 
buffer A and then bound RNA aptamers were eluted from 
the filter using elution buffer B (0.3 M NaOAc, 0.1% 
SDS). The RNA aptamers were reverse transcribed into 
cDNA, amplified by PCR and then RNA was synthesized 
from the synthesized cDNA at the next round of selection. 
The next two rounds of selection comprised purification 
using the MagneHis Protein Purification System 
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(Pr omega). The final round of selection consisted of filtra- 
tion through nitrocellulose filter. The yielded ligand DNA 
was cloned into the pGEM-T vector (Promega) and indi- 
vidual clones were sequenced using the ABI PRISM 310 
genetic Analyzer with the Big-Dye Terminator v3.1 cycle 
sequencing kit (Applied Biosystems). 

Preparation of RNA samples 

Samples of short RNAs (<18nt) were chemically 
synthesized using a DNA/RNA synthesizer (Expedite 
8909, Perseptive). Long RNA samples (>45nt) were 
enzymatically synthesized by in vitro transcription using 
the AmpliScribe™ T7 High Yield Transcription kit. 
Template DNA plasmids for transcription were con- 
structed by inserting the DNA fragment into the 
pGEM-3zf(+) vector (Promega) digested with EcoRI- 
Hindlll and cleaved by Hindlll digestion to enable 
run-off transcription with T7 RNA polymerase. The 
DNA fragments were prepared by annealing primer 
pairs (see the primer list in Supplementary Table SI). 
RNA samples were purified by resolution on PAGE 
under denaturing conditions with 8M urea and 
concentrated by ethanol precipitation. 

EMSA experiments 

Long and short RNAs were mixed in binding buffer A 
containing ImM DTT, 0.3% ribonucleoside vanadyl 
complex and 5% glycerol and in binding buffer B 
(10 mM Tris-HCl, pH 7.5 and 10 mM NaCl), respectively, 
before EMSA with PAGE. The mixtures were reacted for 
30min at 37 and 4°C for the long and short RNAs, re- 
spectively. Incubation temperatures and buffer conditions 
were optimized to establish appropriate binding condi- 
tions and reproducibility between long or short RNAs 
and BsHfqs. Protein-bound long and short RNAs were 
separated from free RNA by 6% or 8% PAGE using 
acrylamide:bisacrylamide (40:1) in 1 x and 0.5 x TEB 
(Tris-borate-EDTA) buffer, respectively. Non-denaturing 
gel electrophoresis proceeded for 60-70 min at lOOV/cm 
on 100 x 100 mm plates at an ambient temperature of 4°C. 
Bands containing RNA on the gels were visualized by 
staining with ethidium bromide or Toluidinblau O 
(Chroma Gesellschaft Schmidt & Co.). 

Footprinting assay 

We prepared 5'- or 3'- 32 P-labeled RNAs (0.07 or 
0.3pmol) in RNA structure buffer (Ambion) containing 
1 ug of yeast tRNA in the absence or presence of 
0.58 pmol BsHfq-His. The samples were incubated at 
37°C for 30 min before digestion with RNases in a 
reaction volume of 10 ul. Each sample was then adjusted 
with either 2 ul of 0.003 U/ml RNase T 2 (Ambion), 2 ul of 
0.0005 U/ml RNase V x (Ambion) or 2 ul of 5 U/ml RNase 
Si (Ambion) at 37°C for 5 min. The reaction samples were 
quenched by adding 20 ul of inactivation/precipitation 
buffer (Ambion) at — 20° C for 1 h, separated by centrifu- 
gation at 20 000 g at 4°C for 20 min, and then pellets were 
washed with 70% ethanol. Dried pellets were dissolved in 
gel loading buffer II (Ambion) and then the digests were 
separated by 10% PAGE under denaturing conditions 



with 8M Urea, and detected by autoradiography 
(Fujifllm). 

UV cross-link 

Samples of 32 P-RNA and BsHfq-His were incubated at 
37°C for 30 min in binding buffer A containing ImM 
DTT and 5% glycerol. Concentrations of the samples in 
the binding buffer A were adjusted to 50 000 cpm/ul for 
the labeled RNA and 0.8 uM for BsHfq-His. The mixtures 
were exposed for 15 min to five 8W germicidal lamps 
(254 nm wavelength) using UV Stratalinker 1800 
(Stratagene) and then incubated with RNase A at 37°C 
for 10 min. Cross-linked samples were analyzed by 15% 
SDS-PAGE and visualized using the Bio-imaging 
Analyzer System (Fujifllm). 

X-ray crystallography 

GST-tagged protein was expressed and purified for crystal 
structure analysis and BsHfq-RNA complexes were 
crystallized as described previously (32). The crystal struc- 
tures of BsHfq-RNA were resolved at BsHfq:RNA molar 
ratios of 1:1 and 1:2 by molecular replacement with a 
model based on the structure of S. aureus Hfq (1KQ2) 
using the program MOLREP (34) in the CCP4 package 
(35) and refined with the program REFMAC in the CCP4 
package and CNS 1.2 (36). Atomic models were fitted into 
electron density maps using the graphics program 
XtalView/Xfit (37) and COOT (38). The quality of both 
crystal structures, assessed using PROCHECK (39), were 
as predicted or better for a structure at this resolution. 



RESULTS 

In vitro selection of RNA bound to BsHfq 

To define a specific RNA sequence that BsHfq recognizes, 
RNA aptamers were obtained by nine SELEX cycles using 
nitrocellulose filter and magnetic beads. Table 1 summar- 
izes the concentrations of BsHfq and RNAs for each cycle. 
Sequencing the 47 clones isolated from the RNA pool 
after nine selection cycles showed that the 22 RNA 
aptamers possessed AG repeats (AG)„ where n > 2. The 
other 25 RNA aptamers did not contain AG repeats. 

To elucidate how BsHfq recognizes AG repeats in RNA 
aptamers, we prepared 25 RNA aptamers with or without 
AG repeats by in vitro transcription using T7 RNA poly- 
merase and investigated their ability to bind BsHfq using 
an EMSA (Supplementary Figure SI). Bands correspond- 
ing to RNA-protein complexes of aptamers containing 
AG repeats shifted, except for aptamer m36f that con- 
tained a (AG) 2 A sequence but did not bind BsHfq. In 
contrast, bands for aptamers without AG repeats did 
not shift (m4f, ml8f and m33f in Supplementary 
Figure SI). These findings indicated that the AG repeat 
plays a key role in BsHfq binding. Figure 1 shows the 
sequences of RNA aptamers with high binding affinity 
for BsHfq in EMSA. Analysis of the aligned sequences 
indicated that all A residues in the sequence of the AG 
repeat were strictly conserved, whereas G residues could 
be replaced with A or U residues. To reconfirm the 
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relationship between AG repeats and BsHfq binding 
ability using EMSA, we used aptamer m49f of which the 
binding affinity of the (AG) 4 A sequence was the highest 
among the aptamers and 28f isolated from the SELEX 
round 6 pool that contained no AG repeats (Figure 2A). 
The amount of the m49f RNA-protein complex obviously 
increased with increasing amounts of BsHfq (Figure 2B). 
The band corresponding to free RNA (m49f) disappeared 



Table 1. Summary of in vitro selection of RNA aptamers 



Cycle RNA (uM) BsHfq (uM) tRNA (uM) 



1 5.0 0.17 

2 2.5 0.08 

3 2.0 0.03 

4 2.0 0.03 

5 2.0 0.02 10 

6 2.0 0.02 20 

7 2.5 0.33 

8 2.0 0.33 20 

9 2.0 0.02 20 



m4 9 f UAGAGAGAGAUUAGAUCCUGUCCGCGGCUA- 

ml 2 f UAGAGAGAAAUUAGUCUAUUGUUGAGCGGC - 

m2 8 f UAGAUAGAGAGAGUCCGAUUUGAACUCGGC 

m3 2f UAGAGAAAGAUUGUCGAAGUUUGUGCACGAU 

m4 0f UGCAGAGAGAGAUAUCCCGGGCCGCGCG 

m4 8f UAAAGAGAGAGAAGUCCGCGGUCGGUCAGU 

m9 f UAGAGAGAGAUAGAUC CGCGGUUUG C CG 

m57f UUGG C C G AGAGAUAGAG AUUUAG C G AAUG C 

m7 f GCUAAUGUGAGAUAUAGAGAGAUAAAAUGU 

ml 4 f UUGAGUAAAGAGAGAGAUGCGAUAUCGCCU 

mlf UGUAGAGAGAGAUAUUCCCGGAUCGGCCAG 

m43f UUG G C C GAG AG AUAG AG AUUUAG C G AUUG C 

m2 0f UGUAUAGAGAUAGAUUC C C AAUGGUC AG C C 

m5f UGAGAGAGAGAGAAACGUUCGUCGACCGUA 



Figure 1. Random sequence regions of selected aptamers with high 
binding affinity for BsHfq determined by EMSA. Names of individual 
clones are left of sequences. Shared sequence motif (AG repeats) is 
shown in bold. Multiple sequences were aligned using ClustalW and 
manually edited. 



at a protein:RNA molar ratio of 1:1 (Figure 2B, lane 10), 
which agreed with the experimental finding that bacterial 
Hfqs form hexamers (2,4,25,30,40). The apparent dissoci- 
ation constant for the complex of m49f and BsHfq 
estimated by the ImageJ program (National Institutes of 
Health, Bethesda, MD) was approximately lOnM. Gel 
bands corresponding to free RNA in the case of 28f 
were not shifted despite the presence of an equimolar 
amount of BsHfq (Figure 2C), indicating the absence of 
affinity for binding BsHfq. These results indicate that an 
AG repeat sequence is necessary for the specific binding of 
BsHfq to its target RNA. 

Enzymatic footprinting assays 

The results of the SELEX experiment suggested (AG)„ 
(n>2) as a candidate sequence for BsHfq binding. To 
identify protein interaction sites in RNA aptamers, 
RNA secondary structure models for several RNA 
aptamers were predicted using the MFOLD program 
package (41). All predicted secondary structure models 
indicated that the AG repeat sequence adopts a loop 
structure (Supplementary Figure S2). We then performed 
a ribonuclease probing experiment using 32 P-5 / - or 3'- 
m49f digested with RNase T 1? RNase Vi and RNase S x 
(Figure 3). The RNase digestion data showed that the AG 
repeat region forms a single-stranded structure 
(Figure 3A, lanes 4 and 8). Figure 3B shows the most 
probable secondary structure of m49f resulting from the 
enzymatic cleavage determined by probing analysis and 
secondary structure prediction. Footprinting analysis of 
32 P-5 / -m49f with BsHfq using RNase T x and Si located 
protected residues in a region containing the AG repeat 
sequence [Figure 3 A (lanes 5 and 9) and B]. Residues in 
this region were also protected by BsHfq in the experiment 
using 32 P-3 / -m49f (data not shown). Consequently, BsHfq 
binds to the loop structure formed by the AG repeat 
sequence. 



A m49f 5' gggacacaauggacgUAGAGAGAGAUUAGAUCCUGUCCGCGGCUAuaacggccgacaugagag 3' 



28f 5' gggacacaauggacgAGGGAUUUGGGUGAUAUUUUUAGUCGGUCAuaacggccgacaugagag 3' 




123456789 10 11 12 1234567891011 12 



Figure 2. Affinity and specificity of BsHfq-binding RNA. (A) RNA sequences isolated by SELEX: m49f (round 9) and 28f (round 6). Nucleotides 
comprising AG repeats highlighted in bold. Nucleotides in lower-case letters indicate constant flanking sequences. EMSA of m49f (B) and 28f (C) 
(1 |iM each) titrated with increasing amounts of BsHfq-His to determine binding to BsHfq. 



1860 Nucleic Acids Research, 2012, Vol. 40, No. 4 



B 



60- 



50 - 



40 - 




(BsHfq 
-His) 



CD 
=3 



mm mm 



A c ^ 



20 

^ V 

G /<b 
G 



\ U 



A A 

/ t 



/ 



G-a-G- 



A 

l<h= 
C 

u \ 



10 - 



123456789 



Ti 




Strong 




Weak 


Vi 




Strong 




Weak 


Si 


^> 


Strong 




=> 


Weak 



Figure 3. Identification of BsHfq-binding site on RNA aptamer m49f. (A) Enzymatic footprinting of BsHfq bound to m49f. Cleavage products of 
32 P-5 ; - m49f modified by RNase T l5 Vi and Si with or without BsHfq-His were analyzed by denaturing gel electrophoresis. Lanes C, AH and G 
represent untreated RNA, alkaline hydrolysis ladder and G ladder, respectively. The vertical line on the right of the gel indicates protected residues 
that are resistant to RNase digestion in the presence of BsHfq-His. (B) Summarized results of enzymatic probing analysis of RNA aptamer m49f. 
Secondary structure was obtained by combining experimental results with MFOLD predictions (41). Major and minor cuts are indicated as symbols 
(see box). Residues corresponding to vertical line in (A) are highlighted in bold. 



Binding affinity of mutant RNAs for BsHfq 

We constructed mutant RNAs based on the sequence of 
aptamer m49f (hereafter referred to as 49R, Figure 4A) to 
determine the importance of residues in AG repeat se- 
quences for BsHfq binding. We initially examined the 
BsHfq-binding affinities of three mutants in which GA 
residues were replaced with CC residues (named 
49R-mtl, 49R-mt2 and 49R-mt3; Figure 4A). Mutant 
RNAs were labeled with 32 P, exposed to short-wavelength 
UV light in the presence of BsHfq and then analyzed by 
SDS-PAGE (Figure 4B). The intensity of bands corres- 
ponding to three 32 P-labeled mutants (49R-mtl, 49R-mt2 
and 49R-mt3) was apparently reduced, compared 
with that of 49 R that contains an (AG) 4 A sequence 
(Figure 4B, lanes 1-4). The EMSA results showed that 
the mutants cannot form stable complexes with BsHfq 
even in the case of 49R-mtl (data not shown). These 
findings indicated that reducing the number of AG 
residue repeats decreased their protein binding affinity. 

The sequence of aptamers obtained from the SELEX 
experiment also indicated that G residues in the AG 
repeat sequence can be replaced by either A or U 



residues (Figure 1). The G residue at the sixth position 
of the AG repeat sequence is in fact, frequently substituted 
in this manner. To determine the nucleotide sequence spe- 
cificity for the binding of these aptamers by BsHfq, we 
generated four more mutants based on the 49 R 
sequence (bottom of Figure 4A); the G residue at the 
sixth position of the AG repeat sequence of 49 R was 
replaced with U, A, C or GG residues and the mutants 
were named 49R-U, 49R-A, 49R-C and 49R-GG, respect- 
ively. The results of UV cross-linking indicated that 
mutant 49R-A and 49 R formed similarly stable complexes 
(Figure 4B, lanes 5 and 7). However, the 32 P signal inten- 
sity of 49R-U, 49R-C and 49R-GG was reduced, 
indicating that these mutants do not easily form 
complexes (Figure 4B, lanes 6, 8 and 9). The results with 
49 R and mutant proteins in EMSA were the same (data 
not shown). Thus, dinucleotide repeat motifs (AG) play an 
important role in the efficiency of RNA aptamer binding 
to BsHfq. 

We constructed a series of RNA fragments by deleting 
residues from the 3 ; -end of the 9-mer 5 / -(AG) 4 A-3 / and 
examined interactions between the fragments and BsHfq 
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using EMSA to analyze crystal structures. The protein 
binding ability of fragments that were shorter than 
7-mer (AG) 3 A, was reduced compared with that of 7-, 8- 
and 9-mer fragments (data not shown). Therefore, the 
minimal RNA fragment suitable for structure analysis 
was defined as a 7-mer (AG) 3 A (hereafter referred to as 
AGr). In practice, EMSA with the AGr fragment showed 
that the intensity of the shifted bands increased with 
increasing protein concentrations (Figure 5, lanes 2-6). 
Furthermore, adding equimolar or excess amounts of 
BsHfq led to significantly weakened bands for free 
RNA. The apparent dissociation constant for the 
binding of AGr to BsHfq estimated by the ImageJ 
software was ~1 uM. The same sequence (AGAGAGA) 
was present in both AGr and 49R-mtl, although 49R-mtl 

A 

49R 5 ' - GGGCGAA [/[7C UAGAGAGAG AU UAGAUC C UGUC C GC GGCUAAA. GC U- 3 ' 

49R-mt1 5 ' - GGGCGAA[/[7C UAGAGAGACCU UAGAUCCUGUCCGCGGCUAAAGC[7- 3 ' 

49R-mt2 5 ' - GGGCGAA[/[7C UAGAGACCCCU UAGAUCCUGUCCGCGGCUAAAGC[7- 3 ' 

49R-mt3 5 ' - GGGCGAA[/[/C UAGACCCCCCU UAGAUCCUGUCCGCGGCUAAAGC[7- 3 ' 

49R-U 5 ' - GGGCGAA [/[7C UAGAGAU AG AU UAGAUC C UGUC C GC GGCUAAA. GC U- 3 ' 

49R-A 5 ' - GGGCGAA [/[7C UAGAGAAAG AU UAGAUC C UGUC C GC GGCU AAA GC U- 3 ' 

49R-C 5 ' - GGGCGAA [/[7C UAGAGAC AG AU UAGAUC C UGUC C GC GGCU AAA GC U- 3 ' 

49R-GG 5 ' - GGGCGAA [7[7C UAGAGAGGAGAU UAGAUCCUGUCCGCGGCUAAAGC[7- 3 ' 

B ^ N # ^ o ^ o eP 

1234 56789 

Figure 4. Effects of RNA mutations on RNA binding to BsHfq. 
(A) Mutant RNA sequences (bold type) constructed for UV 
cross-linking analysis. Tag sequence was added to 3'- and 5'-ends of 
RNA (italic) for in vitro transcription. AG repeat sequence is 
underlined. (B) UV cross-linking analysis of mutant RNA binding to 
BsHfq. 32 P-RNA (50 000 cpm) was incubated with 8.3pmol of 
BsHfq-His and then loaded on gels. Lanes 1 and 5, BsHfq-His and 
49R; lanes 2-4, BsHfq-His and mutants with GA(s) in AG repeat 
sequence replaced by CC(s) (49R-mtl, 49R-mt2 and 49R-mt3); lanes 
6-9, BsHfq-His with mutants in which G in AG repeat sequence was 
replaced by U/A/C/GG (49R-U, 49R-A, 49R-C and 49R-GG). 
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contains an additional sequence outside the AG repeat 
sequence. However, 49R-mtl had lower BsHfq-binding 
affinity than AGr (data not shown). This finding 
suggests that the additional sequence in 49R-mtl inhibits 
the binding of 49R-mtl to BsHfq. 

Escherichia coli Hfq and SaHfq reportedly bind to A 18 
and AU 5 G, respectively (25,31,42). We performed EMSAs 
to determine whether BsHfq can bind to these short RNAs 
(Figure 5). The band corresponding to free RNA dis- 
appeared when the concentration of BsHfq was 4-fold 
higher than that of A 18 (Figure 5, lane 12). This indicates 
that BsHfq can interact with A 18 ; however, the dissoci- 
ation constant for the A 18 binding to BsHfq was 
MOuM, indicating that BsHfq has weaker binding 
affinity for A 18 than for AGr. Moreover, a dissociation 
constant of lOnM has been established between EcHfq 
and Ais (42), indicating stronger binding than that 
between BsHfq and A 18 . In contrast, EMSA showed 
that AU 5 G did not bind BsHfq (Figure 5, lanes 14-18). 
Taken together, these findings demonstrated that the AGr 
fragment possesses significant BsHfq-binding affinity. 



Structure of BsHfq-RNA complex 

We experimentally verified that AGr has higher binding 
affinity than A 18 for BsHfq, as described above. To 
analyze the molecular recognition pattern of BsHfq for 
AG repeats, we determined the structures of the BsHfq- 
AGr complex by X-ray crystallography. Crystals of 
BsHfq-AGr were obtained under molar ratios of 1:1 
(co-crystal 1) and 1:2 (co-crystal 2) (32). The space 
groups of these crystals were 1422 and ¥222 for co-crystals 
1 and 2, respectively. The asymmetric unit of co-crystal 1 
contained six monomers of BsHfq and one molecule of 
RNA, that of co-crystal 2 contained three monomers of 
BsHfq and one molecule of RNA (32). Crystal structures 
of co-crystals 1 and 2 were determined at 2.2 A resolution 
by molecular replacement using a polyalanine model of 
SaHfq (PDB code: 1KQ2) (25), and refined to R/R frQQ 
values of 21.4/22.8% and 25.9/26.0%, respectively 
(Table 2). 



A 18 AU 5 G 

0.05 0.1 0.15 0.2 0.4 - 0.05 0.1 0.15 0.2 0.4 

S I 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 
Figure 5. Affinity of BsHfq for RNAs with sequences AGr, A 18 and AU 5 G. RNA (0.1 mM) was titrated with increasing amounts of BsHfq. 
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Table 2. Diffraction data and refinement statistics for BsHfq-AGr complexes 



Data collection statistics 



Co-crystal 1 



Co-crystal 2 



Date collection 
Space group 

Unit cell parameters (A) 
Wavelength (A) 
Resolution (A) 
Multiplicity 

Number of observations 
Number of unique reflections 
Completeness (%) 



Da" 1 ) 



Average // a(7) 
Matthews coefficient (A 3 
Solvent content (%) 
Refinement 
Data range (A) 
^-factor (%) 

^free (%) 

Number of protein atoms 
Number of RNA atoms 
Number of solvent atoms 
RMSD of bond lengths (A) 
RMSD of bond angles (°) 
Ramachandran Plot 
Most favored regions 
Additional allowed regions 
Generously allowed regions 
Disallowed regions 



1422 

a = b = 123.70, c = 119.13 
1.0000 

25.00-2.20 (2.28-2.20) 

14.8 (14.9) 

350287 

23 732 

100 (100) 

6.3 (32.5) 

38.02 (7.07) 

2.34 

48 

25.00-2.20 

21.4 

22.8 

3333 

154 

138 

0.022 

2.004 

93.1 
6.6 
0.3 
0.0 



F222 

a = 91.92, b = 92.50, c 
1.0000 

30.00-2.20 (2.28-2.20) 

7.4 (7.4) 

92 082 

12 626 

100 (100) 

3.3 (30.7) 

45.75 (4.84) 

2.51 

51 

30.00-2.20 

25.9 

26.0 

1632 

132 

25 

0.035 
2.957 

88.7 
10.2 
1.1 
0.0 



114.92 



£j | /i(hkl) — </(hkl)> | / S hkl Ej/^hkl). Values in parentheses are for the highest resolution shell. 



The quaternary structure of BsHfq in these complexes 
was a homohexameric ring and each subunit of 
homohexameric BsHfq contained an Sm-like domain 
fold consisting of an N-terminal oc-helix followed by a 
five-stranded antiparallel (3-barrel (Figure 6A, B and C), 
as found in other bacterial and archaeal Hfqs (25-30). The 
superposition of these monomer structures between 
co-crystals 1 and 2 yielded a low RMSD value of 0.27 A 
for backbone (N, Coc, C) atoms of secondary structure 
elements. The AGr bound to the distal site of BsHfq 
with a circular conformation (Figure 6), although 
electron densities in a pocket surrounded by F41 and its 
neighboring residues at the proximal site could not be in- 
terpreted. A full-length AGr (7-mer) was identified in 
co-crystal 1. On the other hand, a 3'-end A residue of 
AGr was not found in co-crystal 2; instead, a 5'-end A 
residue from neighboring AGr molecule was found in that 
location. In this manner, we found that although the 
observed length of AGr differed between co-crystals 1 
(7-mer) and 2 (6-mer), the nucleotide binding patterns 
were the same for both crystals except for the 3'-end A 
residue. All riboses in AGr exhibited the C2 f -endo form. 
The A and G residues in AGr adopted anti and syn glyco- 
side bond conformations, respectively (Figures 6A 
and B), although one G residue in co-crystal 2 adopted 
the anti conformation because of intermodular crystal 
contact with Y20 from a neighboring molecule. Both 
structures revealed that A bases in the anti form were 
inserted into the binding pockets formed by two 
aromatic rings of F24 and F29, and that BsHfq was 



recognized by n-n stacking interactions and four 
hydrogen bonding interactions (A Nl - S60 yOH, A N6 
(NH 2 ) - T61 yO, A OH2' - G28 CO and A OH2' - Q30 
sOl) (Figure 6D, E and F). All three G bases in the syn 
form in co-crystal 1 formed a stacking interaction with the 
side chain of Q30 and two of them were involved in a 
hydrogen bonding interaction between 06 of G base and 
rjNH 2 of R32 (Figure 6E and F). The G base in the anti 
form in co-crystal 2 mentioned above displayed neither the 
stacking interaction with Q30 nor the hydrogen bonding 
interaction with R32. The remaining two G bases in the 
syn form in co-crystal 2 revealed stacking interactions with 
the side chain of Q30, but only one of them was engaged in 
hydrogen bonding interaction with the side chain of R32. 
Notably, the atomic coordinates of the side chain of R32 
were not clearly determined due to an ambiguous electron 
density. As a result, except for the hydrogen bonding 
interaction of the G base, co-crystals 1 and 2 represented 
a common interaction mode between dinucleotide repeat 
motifs (AG) and BsHfq, namely two stacking and four 
hydrogen bonding interactions for the A base and one 
stacking interaction for the G base. 

Effects of point mutations in BsHfq on binding to the 
RNA aptamer 

The crystal structures of BsHfq-AGr showed that AGr 
binds to the distal site of BsHfq. On the other hand, the 
electron density was ambiguous in the proximal site of 
BsHfq, rendering feasible the notion that the AG repeat 
sequence also binds to the proximal site. Therefore, 
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Figure 6. Overall structures of BsHfq and the RNA aptamer, AGr. (A) Quaternary structure of BsHfq-AGr in co-crystal 1. One RNA molecule 
binds to BsHfq hexamer. Water molecules are shown as magenta spheres. (B) Quaternary structure of BsHfq-AGr in co-crystal 2. Symmetry-related 
molecule in BsHfq-AGr is shown in light gray and light green. Two RNA molecules bind to one BsHfq hexamer. Water molecules are shown as 
magenta and gray spheres. (C) Superposition of Ca trace of subunits from several bacterial Hfqs. Co-crystals 1 and 2 are shown in black and gray, 
respectively. SaHfq (PDB ID: 1KQ1) is shown in magenta, SaHfq-AU 5 G (PDB ID: 1KQ2) in red, EcHfq (PDB ID: 1HK9) in cyan, EcHfq-poly(A) 
(PDB ID: 3GIB) in blue and P. aeruginosa Hfq (PDB ID: 1U1S) in orange. Pairwise root-mean-square deviations of backbone (N, Ca, C) atoms for 
secondary structure elements in these subunits from co-crystal 1 structure were <0.36A. (D) BsHfq and AGr are shown as gray surface with cyan 
(A residue) or magenta (G residue) sticks, respectively. Part of 2F Q -F C electron density map (yellow) contoured at 1.5 a is shown for AGr. 
(E) Close-up view of BsHfq-AGr interface. BsHfq is shown as semi-transparent molecular surface (light gray). Polypeptide backbone is drawn as 
line representation (light green), while AGr is shown as ball-and-stick representations and colored as in (D). Side chains involved in RNA-binding are 
labeled and are shown as gray bonds colored according to atom type. Hydrogen bonding and stacking interactions between RNA and Hfq are shown 
as dashed lines and arrows, respectively. (F) Schematic representations of interactions between BsHfq and A/G residues in AGr. Dashed lines and 
arrows indicate hydrogen bonding and stacking interactions, respectively. 



we examined whether AG repeats (49 R) specifically 
interact with the distal or proximal site of BsHfq. Based 
on the information gained from RNA-protein interaction 
networks in the distal site of BsHfq and in the proximal 
site of SaHfq (25), we constructed mutant BsHfqs contain- 
ing a single amino acid substitution in the distal (F24A, 
F29A, Q30A and R32A) and proximal (N40A, F41A, 
K56A and H57A) sites (Figure 7A). We then analyzed 
all of the mutant proteins and 49 R RNA using EMS A 
(Figure 7B). Shifted bands corresponding to RNA- 
protein complexes were not evident for mutant proteins 
with substituted residues on the distal site (F24A, F29A, 
Q30A and R32A) (Figure 7B, lanes 3-6), indicating that 
the four mutant proteins lose the ability to bind 49 R 
RNA. This observation suggests that hydrogen bonding 
interaction between atoms 06 of G base and r|NH 2 of 
R32 is needed for AGr recognition by BsHfq (Figure 6E 
and F). Notably, the four residues are highly conserved 
among Gram-positive bacteria (Supplementary 
Figure S3). On the other hand, shifted bands for mutant 



proteins N40A and F41A in which residues on the 
proximal site were substituted, were similar to those of 
the wild-type (Figure 7B, lanes 2, 7 and 8). 
Electrophoretic bands of free RNA were smeared due to 
the dissociation of complexes with mutant proteins K56A 
and H57A in which residues on the proximal site were 
substituted (Figure 7B, lanes 1, 9 and 10). Gel filtration 
chromatography and SDS-PAGE analysis revealed that 
the ability of K56A and H57A to form hexameric struc- 
ture was reduced to ~60 and 30% of that of the wild-type, 
respectively, and that mutant proteins were more sensitive 
to denaturation by SDS than the wild-type (peak $ in 
Figure 7C and D, lanes 2, 4 and 5). Consequently, the 
dissociation of complexes with K56A and H57A was 
caused by disruption of the hexameric formation of 
BsHfq. On the other hand, wild-type and mutant 
proteins, except for K56A and H57A, adopted a 
hexameric structure in solution since their ability to 
form the hexameric structure was maintained and the 
proteins were insensitive to denaturation by SDS 
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Figure 7. Point mutational analysis of distal and proximal sites of BsHfq. (A) Surface representations of distal (left panel) or proximal (right panel) 
sites of BsHfq. Mutated residues are labeled and their side chains are shown in cyan. These figures were generated using PyMOL software. 
(B) EMSA of wild-type and mutant BsHfq binding to RNA (49R). RNA (1 uM) was incubated in the absence (lane 1) or presence of 1.5 uM 
wild-type (WT, lane 2) or mutant (lanes 3-10) BsHfqs. (C) Gel filtration profiles of wild-type (black), F29A (magenta), K56A (cyan) and H57A 
(green) mutants. Peaks corresponding to BsHfq and lysozyme are indicated as $ and #, respectively. (D) A coomassie-stained SDS acrylamide gel 
shows wild-type and mutant BsHfqs. Samples collected from peak $ of gel filtration were resolved by SDS-PAGE. Plus and minus symbols on the 
top of the gel indicate samples boiled or not in SDS buffer for lOmin at 95°C, respectively. Apparent molecular weights of monomeric and 
hexameric BsHfqs are ~9.6 and 58kDa, respectively. 



(peak $ in Figure 7C and D, lanes 2 and 3). Taken 
together, only the distal site of BsHfq possesses an AG 
repeat-containing RNA-binding site. 

DISCUSSION 

Mutagenesis studies of EcHfq and structure analysis of 
the EcHfq-poly(A) complex have revealed that the distal 
site of EcHfq is used to bind to poly(A) (30,31,42). The 
superposition of the hexameric Hfqs of EcHfq-poly(A) 
and BsHfq-AGr results in the RMSD of 0.44 A for 
backbone (N, Coc, C) atoms of secondary structure 
elements. Moreover, studies of the electrostatic surface 
potentials of EcHfq (30) and BsHfq have revealed that 
the distal and proximal sites of these Hfqs are all positive- 
ly charged (Supplementary Figure S4). However, the se- 
quences and structures differ between RNA-binding sites 
at the distal site of EcHfq-poly(A) and BsHfq-AGr, and 
they produce distinct binding patterns, as shown below. 
Firstly, Hfq-bound RNAs in the BsHfq-AGr confor- 
mationally differ from those in EcHfq-poly(A). At least 
both the R and A sites (where R is a purine nucleotide) are 
necessary for poly(A) recognition in EcHfq (30), and a 
residue at an N site (where N is any nucleotide) is there- 
fore likely to be essential to connect to nearby residues at 
R and A sites. Hence, EcHfq might recognize A-R-N 
trinucleotide as the smallest unit of poly(A) (right panel 
in Figure 8A). On the other hand, we demonstrated that 



residues at A (F24 and F29) and G (Q30 and R32) sites are 
absolutely required for specific RNA binding at the distal 
site in BsHfq (Figure 7). Thus, BsHfq interacts with AG 
dinucleotide as the smallest unit of the AG repeat motif 
(left panel in Figure 8A). As a consequence, both RNA 
shapes adopted a circular conformation, whereas poly(A) 
in EcHfq assumed a waving conformation, and then the N 
base (the base in the N site) in EcHfq is exposed to the 
solvent (Figure 8). The RNA-protein binding patterns 
differed between BsHfq-AGr and EcHfq-poly(A). The 
G base at the G site of BsHfq formed a hydrogen 
bonding interaction with the side chain of R32 (left 
panel in Figure 8B). Mutational analysis also revealed 
that R32 is critical for interaction with the AG repeat 
motif (Figure 7). The residue corresponding to R32 is 
Q33 in EcHfq (right panel in Figure 8B), and the side 
chain of Q33 does not interact with the N base in EcHfq 
(30). Instead, the interaction of Q33 of EcHfq with the A 
base (the base in the A site) is important for poly(A) 
binding (30), although this binding pattern does not 
exist in BsHfq (Figure 8B). Notably, R32 of BsHfq and 
Q33 of EcHfq are highly conserved among Gram-positive 
and Gram-negative bacteria, respectively (Supplementary 
Figure S3). These findings indicate that the binding 
pattern of residues at G site to target RNA is specific to 
BsHfq and not EcHfq. Lastly, both the A base of BsHfq 
and the R base (the base in a R site) of EcHfq are inserted 
into the pocket formed by residues corresponding to F24 
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Figure 8. Molecular recognition of RNA by distal site of Hfq. (A) Left panel: surface representation of hexameric BsHfq (gray) with stick repre- 
sentation of A (cyan) and G (magenta) residues of AGr. Right panel: surface representation of hexameric EcHfq (gray) with stick representation of 
poly (A) (cyan). Subunit boundaries are indicated by dashed lines. (B) Stereo view of BsHfq-AGr and EcHfq-poly(A) (left and right panel, 
respectively). Molecules are shown using same representation and color scheme as shown in Figure 6E. R and N indicate purine nucleotide and 
any nucleotide, respectively. These figures were generated using PyMOL software. 



and F29 of BsHfq. However, distinct binding pattern of 
hydrogen bonding interactions in BsHfq and EcHfq 
directly influence specific structural features inside the 
pocket. The A base of BsHfq forms hydrogen bonding 
interactions with S60 and T61 (Figures 6E, and F and 
8B), whereas the R base of EcHfq does so with Q52 and 
T61 (30). In particular, the S60 residue of BsHfq locates 
near the bottom of the pocket, whereas the Q52 residue 
locates outside the pocket (Figure 8B). Thus, only the 
binding pattern of S60 in BsHfq allows insertion of the 
A base deeply into the pocket without steric clashes. These 
results show that the distal site of BsHfq can selectively 
bind to AGr in a sequence- and/or structure- specific 
manner. 

Co-crystals 1 and 2 had some ambiguous, but 
non-protein electron density at the proximal site, especial- 
ly in a pocket surrounded by F41 and its neighboring 
residues (Supplementary Figure S5). The electron densities 
cannot be interpreted as AGr, because our mutational 
analysis revealed that AGr does not interact with the 
proximal site of BsHfq (Figure 7). Moreover, we con- 
firmed that the BsHfq used in this study contained no 
residual nucleic acid after purification by the Warburg- 
Christian method (43). A positively charged region in 
the proximal site of SaHfq appears to interact with 
small U-rich RNA (AU 5 G) (25) (Supplementary 
Figure S4). Moreover, AU 5 G-bound amino acid residues 
in the proximal site were identified by the structure of 
SaHfq-AU 5 G complex (25). Structure-based sequence 



alignment revealed that AU 5 G-bound amino acid 
residues and the conformations of their side chains in 
SaHfq are also conserved in BsHfq (Supplementary 
Figures S3 and S5). Thus, BsHfq is likely to bind to 
AU 5 G, but in fact, AU 5 G did not bind to BsHfq 
(Figure 5). On the other hand, the crystal structure of 
Salmonella typhimurium Hfq (StHfq) in complex with U 6 
RNA has recently been determined, in which U 6 RNA 
interacts with the proximal site of Hfq (44), and a free 3' 
hydroxyl group of RNA is crucial for high-affinity binding 
to StHfq. In particular, the S'-terminal guanosine in RNA 
obviously reduced StHfq binding. Thus, we speculate that 
the absence of interactions between AU 5 G and BsHfq 
could be due to the guanosine residue at the 3'-end of 
AU 5 G. Therefore, these findings raise the notion that 
BsHfq still associates with some RNA(s) through the posi- 
tively charged proximal site, as found in SaHfq or EcHfq 
(25,31) (Supplementary Figure S4). However, whether or 
how the proximal site of BsHfq recognizes RNA(s) 
remains unclear. 

Hfq in Gram-negative bacteria functions as a 
post-transcriptional regulator that acts by mediating inter- 
actions between many sRNAs and their mRNA targets. 
On the other hand, BsHfq and SaHfq that are both from 
Gram-positive bacteria interacted with different RNA 
sequences (mRNA and/or sRNA) via structurally 
distinct RNA-binding surfaces (distal or proximal site). 
BsHfq-binding sequences in SRI sRNA and ahrC 
mRNA have been elucidated (24), and they are involved 
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in short AG repeats GAAUAAGAGA (SRI sRNA) or A 
AAUAGAG (ahrC mRNA). Therefore, these results when 
combined with the findings of the SELEX experiments 
indicated that the target sequence for BsHfq certainly 
comprises AG repeats. Taken together with these data, 
the results presented herein suggest that SRI and ahrC 
mRNA recognize the same distal site of BsHfq. We 
further examined whether the intergenic regions of 
B. subtilis 168 genome include the AGr sequence using 
the SubtiList database (http://genolist.pasteur.fr/ 
SubtiList/). We identified the sequence in 110 intergenic 
regions and in 215 protein coding regions, indicating that 
the AGr sequence does in fact exists in the B. subtilis 
genome. 

Amino acid residues on the surface of the distal 
site, including those that bind to AGr in BsHfq, are 
highly conserved in Gram-positive bacteria such as 
S. aureus and L. monocytogenes (Supplementary 
Figure S3). Moreover, lmo0850 mRNA, which binds to 
L. monocytogenes Hfq, also possesses AG repeats (45). 
Therefore, we speculate that several other Hfqs might 
function through interaction with the AG repeat motif 
among Gram-positive bacteria. To date, various RNA 
target sequences have been identified and they can dis- 
criminate the RNA-binding surface(s) of hexameric Hfq 
at the atomic level. Furthermore, multiple functional roles 
of Hfq have been identified in several bacterial species, 
probably as a consequence of this. Hence, bacterial Hfqs 
might have species-specific functions among bacterial 
phyla. To understand Hfq activities, recognition patterns 
between RNA(s) and Hfq across various biological species 
should be further elucidated in future studies from both 
structural and functional perspectives. 
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